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Abstract: 3D Georgi-Glashow model is studied on the lattice in the London limit 
in an infrared but an intermediate region before the screening appears. Abelian and 
instanton dominances are observed after abelian projections in a unitary gauge and 
roughly in the maximally abelian gauge. Using an inverse Monte-Carlo method, we 
determine an effective instanton action in both gauges. When we restrict ourselves to 
some regions of parameters j3 and k, we obtain an almost perfect instanton action, 
performing a block-spin transformation on the dual lattice. It takes a form of a 
Coulomb gas and reproduces fairly well the string tension obtained analytically by 
Polyakov. The almost perfect actions in both gauges look the same in the infrared 
region, which suggests gauge independence. 
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1. Introduction 

It is very important to understand confinement mechanism of QCD. Wilson's lattice 
formulation [1] shows that the confinement is a property of a non-Abelian gauge 
theory of strong interaction and there are a lot of numerical lattice calculations 
showing the confinement of color. But the mechanism of confinement is still not well 
understood. 

One of approaches to the confinement problem is to search for relevant dynamical 
variables and to construct an effective theory in terms of these variables. In 1970's, 
it was pointed out that the confinement of quarks may be explained as the dual 
Meissner effect due to condensation of monopoles [2, 3]. 

In (QCD) 4 , 't Hooft[4] proposed an idea of abelian projection. After a partial 
gauge fixing (called as abelian projection), SU(N) gauge theory can be reduced to an 
abelian U(l) Ar_1 theory with N — 1 different charges and magnetic charges. Namely, 
N — 1 monopoles appear after the abelian projection. Quark confinement could 
be understood by the dual Meissner effect due to condensation of these monopoles. 
Actually there have been many numerical data supporting the above conjecture when 
one performs the abelian projection in the Maximally abelian (MA) gauge. 

The first interesting findings are some phenomena called abelian dominance. 
Abelian Wilson loops composed of abelian link variables alone seem to reproduce 
essential features of quark confinement in the Maximally Abelian gauge [5]. 

The second is monopole dominance. It has been shown that monopoles alone 
are responsible for abelian quantities in the infrared region. Abelian Wilson loops 
and abelian Polyakov loops are written as a product of monopole and photon contri- 
butions. The monopoles which are defined by DeGrand and Toussaint[6] reproduce 
the confinement features [7]. 

Furthermore an effective lattice monopole action was derived by Shiba and 
Suzuki [8] in a simple case with two-point interactions alone. The work was extended 
to include 4 and 6 point interactions in [9] and [10]. 

There is another system having such a magnetic quantity naturally. It is the 
three-dimensional Georgi-Glashow model (GG) 3 where there exist a famous 't Hooft- 
Polyakov instanton [11, 12] having a magnetic charge. Polyakov showed analytically 
[13] that the string tension of 3D Georgi-Glashow model has a finite value. He made 
a quasi-classical calculation using a dilute Coulomb gas approximation of 't Hooft- 
Polyakov instantons. However justification of his method is not clear. If the quasi- 
classical assumption using the dilute instanton gas is true, abelian and instanton 
dominances should hold in the infrared region of three-dimensional Georgi-Glashow 
model. It is expected naturally that an effective abelian instanton action works well 
in the infrared region. 

It is shown in [14] that Polyakov's picture of the dilute Coulomb instanton gas is 
too simple to describe all infrared physics of (GG) 3 . Actually, it can not explain the 
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screening of doubly charged particle due to W bosons. However, such a screening 
appears only for particles with even charge and also in the long infrared range of 
0(2M w /a), where M w is a mass of the W boson and a is the string tension. Hence if 
we restricted ourselves to an infrared but an intermediate region before the screening 
appears, we may be allowed to adopt such a picture as Polyakov's one. 

It is just the purpose of this study to test the Polyakov's assumptions and to 
check calculations done by him quantitatively formulating the theory on the lattice 
and restricting ourselves to the intermediate region before the screening appears. 
Our method is Monte-Carlo numerical simulations without ad-hoc assumptions. In 
Section 2, we briefly review the Polyakov's work in 1977[13]. We calculate the total 
energy of 't Hooft-Polyakov instanton in the London limit. In Section 3, we show the 
method of our numerical study. We derive a lattice instanton action in 3D Georgi- 
Glashow model using an inverse Monte-Carlo method. We then perform a block- 
spin transformation of DeGrand-Toussaint instantons. These discussions are almost 
similar to those in 4D pure QCD in MA gauge fixing[8]. In Section 4, we show our 
numerical results of the lattice monopole action in 3D Georgi-Glashow model. Then 
we calculate the string tension derived by Polyakov using the parameters determined 
in our numerical simulations. Section 5 is devoted to concluding remarks. 



2. Confinement mechanism in 3D Georgi-Glashow model 

We review briefly the 't Hooft-Polyakov instanton and confinement mechanism due 
to the instantons as discussed by Polyakov in Ref[13]. 

2.1 Georgi-Glashow model and 't Hooft-Polyakov instanton 

The three-dimensional Georgi-Glashow model is given by 

S = j d'xS^F^ + i W + - v 2 ) 2 } , (2-1) 

where 

FZ„ = d,Al - d v Al + ge« hc A\A% (2.2) 

= d,<p a + ge^A^. (2.3) 

It is a SU(2) Yang Mills theory with an adjoint Higgs field. It has a phase 
with spontaneously broken symmetry (SU(2) — > U(l)) and the gauge field becomes 
massive because of the Higgs mechanism. The gauge boson mass is given at the tree 
level by 

m w ~ gv. (2.4) 



2 



In the case of three space-time dimensions, the theory has an instanton solution 
- 't Hooft-Polyakov instanton. One instanton solution takes the form [11] 



r = ^H(gvr), (2.5) 
A l = - g Sa^[l-K{gvr)], (2.6) 

where H and K are dimensionless functions determined by the equations of motion. 
They satisfy 

K(gvr) — > 
H(gvr) — > gvr 



r -> oo. (2.7) 



For the ansatz given in Eqs.(2.5) and (2.6), the energy of the one-instanton 
system is given by 
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(2.8) 



where £ = gvr. This can be calculated numerically solving equations of motion for 
K and H. 

The value of S at this minimum is expressed as 

Sci = ^e(\/g). (2.9) 

g 2 

The function e(x) is slowly varying. For example, it is known [11] that e(0) = An. In 
our case, we calculate Eq.(2.8) in the London limit (A — > oo) because we adopt the 
London limit in lattice calculations done later. We obtain 

e(oo) = 4tt x 1.67. (2.10) 

2.2 Instanton dilute gas model and Confinement 

Polyakov showed that, when the Higgs mass is large enough, the infrared structure 
of the model in its Higgs phase is dominated by a dilute Coulomb gas of 't Hooft- 
Polyakov instantons with long-range interactions (r ^> m w ~ l ) which disorder the 
system. 

When we consider quantum fluctuations around the instanton background, we 
need to consider zero mode solutions the number of which are equal to that of 
conservation laws violated by the classical solution [13]. 
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The partition function in the one-instanton case is given in the tree approxima- 
tion by 1 

Z x ~ J dRN 3/2 m 3 w e- Scl , (2.11) 

where N is the normalization factor 

N = tr f d'x{(F*)* + (T>***)*} 



4:7TV r ^ 1 



9 J £ 



2 



= r ^a(W). (2.12) 

It is evaluated in the London limit (A — > oo) as 

a(oo) ~ 4tt x 3.33. (2.13) 

Let us next discuss a system of several instantons. In Ref. [13], Polyakov consid- 
ered a case in which the positions of instantons are far from each other ( r ^> m^ 1 ) 
and the instantons have only the Coulomb interaction. This is justified when the 
Higgs mass is large. Then the action is given by 

m w 2\\- 2 , 27r „ q a q b rt1 
S = —e{\ g )2^q a + — l^T^ 2 - 14 

g 2 a g 2 t^ b \xa-x b \ 

where q a is an instanton charge and q a = ±1 alone are taken into account. 
The partition function of the dilute gas model is given by 

Z = E ^ / rid^expl-^V ,^^, }, (2.15) 

where 

C = N 3 / 2 m w 3 e-^ e(X/92) . (2.16) 

This is rewritten as 

z = l v * e S/ ( n d3 ^ exp (-K£) 2 I " d3 ^x) 2 +^E^a)| 

n,{q a } ' J=l ^ a J 

= |p X exp|-i(|-) 2 |d 3 a; ((^x) 2 -2M 2 cos X )|, (2.17) 



1 Polyakov evaluated the integral beyond the tree level (Ref.[13]). The one-instanton par- 
tition function in the one-loop approximation is given 
5 2 [$ cl , [$ d) ] ]) cxp (- E rl# o log ^)™^- Sd - 
But we focus only on the tree-level result for simplicity. 
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where 



M 2 = 2\^j (. (2.18) 

The static potential between a pair of static electric charges is estimated by the 
Wilson loop (Fig.l) defined by 

W[C\ = ( e i9 $c A l Ax ») = ( e i9 fs B » AS »), (2.19) 
where the magnetic field from an instanton is 

B,{x) = ^. (2.20) 

One easily recognizes in (2.20) the Dirac quantization condition. In the case of several 
instantons, the strength of the magnetic field is given by 



B,(x) = -f d*y^%k(y), (2.21) 
g J \x — y\ 6 

where the instanton density is 

k(x) = Y / q a S 3 (x-x a ). (2.22) 

a 

The Wilson loop Eq.(2.19) is expressed as 



vv = \e J / = 

where 



W = ( e </d s «/W*W) = (2.23) 



Z[r,\ = J Z^exp j-i- d 3 x[(^(x-^)) 2 -2M 2 cos X ]|, (2.24) 

^)=J 8 ^-^w- (2 - 25) 

The path integral (2.24) may be evaluated by the saddle-point solution (see 
Appendix A): 



W[C] ~ exp{-RTx 



(2.26) 
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It shows an area law in the approximation of large T and large R. The string tension 
is obtained explicitly as 

Mg 2 



Ccl = 
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(2.27) 



5 



time A 




s 



T 



space 



R 



Figure 1: Wilson loop 



Finally in this section, we show that Eq.(2.14) is rewritten in the outside of the 
instantons as 



We will compare later this form of action with a lattice instanton action derived 
numerically. 

3. Lattice Georgi-Glashow model and instantons 

Now let us make an analysis of the Georgi-Glashow model in the framework of 
lattice field theory without recourse to the quasi classical approximation. First we 
check the abelian dominance and the instanton dominance in the infrared region 
numerically. Next we try to derive an effective instanton action in the continuum 
limit and compare it with the action derived by Polyakov (2.28). 

3.1 Lattice Georgi-Glashow model 

The lattice form of Eq.(2.1) is expressed by 




(2.28) 




+ kJ2^(®l(x)®l(x) - ®l{x)U^x)<$> l {x + ii)UI{x)) 



(3.1) 



x 



where the plaquette gauge variable is 



U^(x) = U»(x)U v (x + M(x + v)Ut{x). 



(3.2) 
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In the case of three space-time dimensions, the relations between the continuum and 
the lattice quantities are 

U^ax) = ex.p[iagA°(x)—}, ® L (ax) = a0 a (x)y, (3.3) 

m 2 L = a 2 v 2 \, \l = a\, (3.4) 

(3 = -4-, k 2 = v 2 a. (3.5) 
g 2 a 

For simplicity, we consider only the London limit — > oo. Then 

2ti($ L (x)$ L (x)) — > 1 (3.6) 

and 

E^l = 1- (3-7) 

a 

The action is reduced to 

S = P E (i-^tr^(x))+«]^ 

(3.8) 



3.2 The method 

The method of our numerical study is the following: 

1. We generate thermalized vacuum configurations {U^(x)}, adopting, for sim- 
plicity, a unitary gauge condition 

tf, = 5 a3 . (3.9) 

In this gauge, the Higgs field disappears in the action. There remains still U(l) 
symmetry under the unitary gauge condition and hence it is regarded as one 
of abelian projections. 

2. The lattice considered is 48 3 for (3 = 4.0 ~ 6.0 and k = 1.00 ~ 1.30. We also 
used 12 3 , 16 3 , 24 3 and 32 3 lattices to study volume dependence. 

3. We next extract abelian link variables (m m (x)} in the unitary gauge as 

U^x) = C^x)u^x) (3.10) 

= GU*>( o e -.»„wj- ("D 

4. We generate also the vacuum configurations in the maximally abelian (MA) 
gauge from those in the unitary gauge. They are generated after a gauge trans- 
formation in such a way as maximizing a quantity R = J2x,fi tr (o- 3 U fl (x)a 3 Uj l (x 
The abelian link variables in MA gauge are extracted similarly as in (3.11) . 
The Higgs field reappear in MA gauge, but we don't need them, since we 
evaluate here only operators composed of abelian link variables. 
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5. An integer instanton charge is defined from the U(l) plaquette variables fol- 
lowing DeGrand and Toussaint [6]. The U(l) plaquette variables are written 
by 

e^(x) = e^x) + e u ( x +^)- e^x +»)- e v {x), (-4tt < e^(x) < An). (3.12) 

It is decomposed into two terms: 

Qfj,v(x) = Q^v{x) + 2nn^(x), (-7T < Q^ u (x) < vr). (3.13) 

Here, Q^, v {x) is interpreted as the electro- magnetic flux through the plaquette 
and the integer n liV {x) corresponds to the number of Dirac string penetrating 
a plaquette. One can define a quantized instanton charge as 

11 

= 2~7r ' 2 € ^ vp ^^ up ^ 

= —^iwpd^n vp {x), (3.14) 

where d denotes the forward difference on the lattice. We then obtain vacuum 
ensembles {k(x)} from the vacuum configuration {u^x)} both in unitary and 
MA gauges. 

6. We determine an effective instanton action S[k] from the vacuum ensemble 
{k(x)} with the help of an inverse Monte-Carlo method first developed by 
Swendsen[15] and extended to closed monopole currents by Shiba and Suzuki 
[8]. Here instantons are a site variable and the original Swendsen method works. 
We consider a set of independent instanton interactions that are summed up 
over the whole lattice. We denote each operator as Si[k]. Then the instanton 
action may be written as a linear combination of these operators: 

S[h]=Yl G iSi[k}, (3.15) 

i 

where Gi are coupling constants. 

We determine the set of couplings Gi from the instanton ensemble {k(x)}. 
Practically, we have to restrict the number of interaction terms. It is natural 
to assume that instantons that are far apart do not interact strongly and to 
consider only short-ranged interactions of instantons. This assumption is not 
justified when the screening appears. Non-local long-range interactions are 
inevitable if we express the theory in terms of abelian quantities alone. But 
here we are restricted to the intermediate region before the screening appears. 

The form of actions adopted here is 10 quadratic interactions : 

S[k] = Go E k (x? + GiJ2 k (xMx + H) + G 2 £ ... + C 9 - (3-16) 

X X,fl 

= ]T k(x)D(x - x')k(x'). (3.17) 

x,x' 
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Figure 2: Operators of instanton action 



The detailed form of interactions is shown in Fig. 2. 

7. To study the continuum limit, we perform a block-spin transformation in terms 
of the instantons. We adopt n = 1,2,3,4,6,8 extended instantons which are 
defined as 

71-1 

k( n \ x {n) )= k(nx (n) +d+]2 + k3). (3.18) 

i,j,k=0 

This is a block-spin transformation on the dual lattice and the renormalized 
lattice spacing is b = na((3, k). The value of the blocked instanton k^ n \x^) is 
integer. 
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We determine the effective instanton action from the blocked instanton ensem- 
ble {k( n \x^)} using the above inverse Monte-Carlo method. Then one can 
obtain the renormalization flow in the coupling constant space. 

4. Results 

4.1 String tension 

We first measure the static potential between a pair of quarks and the lattice string 
tension for (3 = 4.0 ~ 6.0 and k = 1.00 ~ 1.30. They are measured from T — > oo of 
R x T Wilson loops. The Coulomb potential is not 1/R but log/? because we are 
considering the three-dimensional case. The static potential is fitted as A log R + 
a L R + B (where A, B are constants). The usual SU(2) Wilson loops and the string 
tension are evaluated to fix the scale of this theory. 

Abelian Wilson loops in terms of the abelian link variables 9 p (x) in Eq.(3.11) are 
measured similarly. The string tension is obtained from the abelian Wilson loops. 
They are rewritten as 

W[C] = exp{i£0 M (z)J M (z)} = exp{~e^(x)M^(a;)}, (4.1) 

where J p (x) = d'^M pu (x) . Then the abelian Wilson loop operator can be decomposed 
into instanton and photon parts as shown in Ref.[7]: Using an identity 

M pv {x) = -Al\x - x')[& a {d^M au - d v M ap ) + ^e a ^e Xpa d' a d x M pa ], (4.2) 
we get 

W[C] = W p [C}-W m [C], (4.3) 

W P [C] = expi-t^dft^Ai'ix-x^Ux')}, (4.4) 

x,x' 

W m [C] = exp{2mJ2Hx)Al 1 (x - x')^e apa d a M pa }, (4.5) 

x,x' 

where A~[ 1 (x) is the lattice Coulomb propagator. 

In Table. 1 and in Fig. 3 we show the data of the string tensions for the regions 
of f3 and k. It is interesting that the instanton contribution to the abelian Wilson 
loop gives us a static potential which is almost linear: V(R) ~ a^R + B. The string 
tension is determined also from this potential. Except for the data of instanton 
contributions in MA gauge which are smaller than the SU(2) ones, we see abelian 
and instanton dominances are fairy good. This supports the analysis of Polyakov[13]. 
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Table 1: Lattice String tension of Georgi-Glashow model 
Potential (Unitary gauge) Potential (MA gauge) 
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Figure 3: q-q potential 



4.2 Instanton vacuum in both gauges 

Polyakov[13] assumed the dilute instanton gas approximation. Let us see the in- 
stanton vacuum configurations. As we mention in the previous section, we calculate 
the instanton configuration using the DeGrand-Toussaint definition. The instanton 
configuration and anti-instanton configuration in unitary gauge condition for (3 = 4.0 
and k = 1.30 are shown in Fig. 5. Those in MA gauge condition are shown in Fig. 6. 
These configurations look to be a dilute gas for this set of coupling constants. 

Polyakov assumed that instanton charges q a = ±1 alone contribute to the sum- 
mation. We plot a histogram of numerically obtained instanton charges in Fig. 7 for 
various blocked instantons. Theoretically an n blocked instanton can take a charge 
up to 3n 2 — 1. However as seen from Fig. 7, the number of instantons having a charge 
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Figure 4: Lattice String tension of Georgi-Glashow model 




Figure 5: Instanton (left) and anti-instanton (right) vacuum configuration (unitary gauge) 




Figure 6: Instanton (left) and anti-instanton (right) vacuum configuration (MA gauge) 



larger than one is very small. This is consistent with Polyakov's assumption. 



4.3 Parameter regions of (3 and k 

Let us discuss here what values of the parameters (3 and k should be chosen. We 
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Figure 7: Instanton charge histogram 



plot the static potential in Fig.8 for different values of re with fixed (3 = 5.0 and the 
instanton density in Fig. 9 for f3 = 6.0. 

We see that, in the case of the unitary gauge condition, the instanton density 
of the system is very large for small re, whereas in the case of MA gauge condition 
it remains low. The dilute gas approximation may be wrong for small re regions in 
unitary gauge. Actually, the abelian and the instanton potentials have large errors 
in the unitary gauge, whereas in MA gauge these potentials are seen clearly. The 
Georgi-Glashow model is reduced to pure SU(2) gauge theory in the limit of re — > 0. 
Such density and potential behaviors of both gauges are seen also in four- dimensional 
pure SU(2) gauge theory [16]. 

What happens in the large re region? One can find that there seems to be a 
critical re c above which the instanton density becomes vanishing and the instanton 
contribution to the potential becomes zero. Note that the Georgi-Glashow model in 
three dimensions is always in the confinement phase contrary to the model in four 
dimensions. The existence of re c does not mean that instantons are not responsible 
for confinement. In the limit of re — > oo, the three-dimensional Georgi-Glashow 



13 



model becomes the three-dimensional compact U(l) model. The confinement in the 
latter model is proved to be due to instantons for all f3 region [17]. We can define 
instantons in the compact U(l) model naturally without following DeGrand and 
Toussaint (DT). It is those instantons which are responsible for confinement in the 
compact U(l) model. (We call it as natural instanton.) This concludes that the DT 
definition of instanton is not good for larger k region than k c in three-dimensional 
Georgi-Glashow model. In Ref.[8], the effective actions of DT and natural monopoles 
are studied in 4D compact QED. They are not always in agreement. Hence it is not 
understandable that the difference appears for such large k and (3 region also in 3D 
(GG)3. We admit that the physical meaning of k c is not understood theoretically, 
however. In this numerical study, we have to adopt the DeGrand- Toussaint definition 
of instanton, since we do not know how to define natural monopoles in 3D (GG)3. 

Instantons in the unitary gauge are more interesting, since they correspond to 
the instantons studied by Polyakov[13]. Hence in the following discussions we restrict 
ourselves to the value of k around 1.0 which is near but less than k c . 

What about the (3 dependence? We show in Fig. 10 the (3 dependence for fixed 
k = 1.00. For small (3, the static potential for large distance is not fixed clearly. This 
comes from that the theory in the small f3 region has a large lattice string tension 
which corresponds to a large lattice distance as shown later. To study the continuum 
limit, it is better to study larger (3 region. 

However, for very large (3, even k = 1.00 becomes larger than k c in MA gauge. 
Namely k c depends on the value of (3. When (3 becomes larger, k c becomes smaller. 
k c in MA gauge is a little bit smaller than that in the unitary gauge. Since we discuss 
around k — 1.00, we are restricted to (3 = 4 ~ 7 regions. 

4.4 Lattice instanton action 

Since the instanton dominance is seen numerically, an effective instanton action is 
expected to work in the infrared region. Let us derive the lattice instanton action 
from instanton configurations using Swendsen's method [15] (See also Appendix B). 
We express D(x) in (3.17) as 

D(x - x') = G 5 3 (x - x') + G^S^x - (x' - i)) + ...) + G 2 ..., (4.6) 

where we restrict ourselves to ten different two-point couplings shown in Fig. 2. Our 
numerical results are the following: 

1. The inverse Monte-Carlo method works very well and the 10 coupling constants 
are fixed very beautifully. The convergence is very fast. 

When we start from the ensemble of blocked instanton configurations {k n (x)}, 
we can get also a lattice action S[k n ] very rapidly. 

We have tried to include a four-point self interaction in addition. In the case 
of n — 1, we can fix the coupling constants, although the convergence is very 
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Figure 8: k dependence of g-g potential 



slow. The convergence was not obtained in the case n > 1. Hence we take into 
account only the two-point interactions. 

We show in Fig. 11 an example of coupling constants in both gauges, where the 
horizontal axis stands for the distance between two instantons in each two-point 
interaction. The data points of 5*8 and Sg in Fig. 11 are almost degenerate. 
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Figure 9: Instanton density 
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Figure 10: f3 dependence of q-q potential 

The action (2.28) derived by Polyakov[13] is composed of the self-interaction 
and the Coulomb term. Hence let us check if the above action obtained nu- 
merically can be written similarly as in (2.28). When we rewrite the lattice 
Coulomb propagator as 



Al\x - x') = C 6 3 (x - x') + d^Or - (x' - 1)) + ...) + C 2 ..., (4.7) 
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we get a beautiful fit 

d ~ Const, x d 0) 
in both gauges as shown in Fig. 11. 



(4.8) 



The lattice instanton action of Georgi-Glashow model in three-dimension can 
be written as 



S = Self term + Const, x J2 k(x)Al 1 (x - x')k(x'). 



(4.9) 



Although the Coulomb term is reproduced well, we can not compare the self 
terms of the lattice and the Polyakov actions, since we need to know the details 
of short-ranged renormalization effects. Actually the lattice Coulomb propa- 
gator is finite even at zero distance, whereas the propagator in the continuum 
is infinite. 
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Figure 11: Coupling constant of lattice instanton action and Coulomb propagator 



3. To study the continuum limit, we need to take both V — > oo and a — > limits. 
First we study the volume dependence of the action, adopting 12 3 , 16 3 , 24 3 , 32 3 
and 48 3 lattices. Fig. 12 is an example of the lattice volume dependence of 
couplings for (3 = 4.0, k = 1.225 and n — 2. In both gauge conditions, the 
lattice instanton action is almost independent of the lattice volume. 

4. We next study the a — > limit, performing the block-spin transformation 
(3.18). To get the renormalization flow, we have to fix a scale unit. Since 
there is no physical experimental data in three-dimensional Georgi-Glashow 
model, we fix the scale in terms of the string tension determined from lattice 
simulations. The lattice spacing is defined with 

a = ./^S (4.10) 
V °phys 



17 




Figure 12: Volume dependence of Coupling constants 



We perform the block-spin transformations using the scale. An example of the 
renormalization flow is plotted on the G\(b) — G 2 (b) and the G 2 {b) — G 3 (b) 
projected planes in Fig. 13. We see that there seems to exist a universal flow 
line especially in the infrared large b = na((3, k) region. Each flow line for 
large b is straight with small errors. The universal flow lines of both gauges 
are almost the same for large b. This suggests the gauge independence of the 
action in the infrared regions of the continuum limit. 




Figure 13: Flow of Coupling constants 



5. To check if the universal curve is the renormalized trajectory which is the 
continuum theory, let us study a scaling behavior. The continuum limit of the 
lattice theory is realized when n — > oo and a — ■> for fixed b is taken. Effective 
actions are usually a function of the step n of the block-spin transformation 
and the lattice distance a(f3, k). If the effective action shows a scaling, that is, 
a function of the product b = na(/3, k) alone, it is the renormalized trajectory 
(RT). Since we have seen the couplings except the self one are well reproduced 
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by the Coulomb interaction, let us concentrate on the constant coefficient of 
the Coulomb term in (4.9). 

The theory has two parameters (3 and k on the lattice and g and v in the 
continuum. For the purpose of the comparison with the continuum theory, 
we tune (3 and k so that we may get the same value of v — \J n/a((3, tz) as in 
Table. 2. Here we fix v ~ 2.56. The values of the parameters (3 and k in Table 
1 are chosen similarly. 



p 


K 


V 


K 


V 


K 


V 


6.0 


1.030 


2.549(68) 










5.8 


1.030 


2.552(61) 


1.005 


2.419(62) 


1.015 


2.519(65) 


5.6 


1.035 


2.554(66) 










5.4 


1.055 


2.566(62) 


1.050 


2.443(67) 


1.065 


2.656(65) 


5.2 


1.065 


2.567(63) 


1.065 


2.468(68) 


1.080 


2.588(58) 


5.0 


1.065 


2.549(66) 


1.095 


2.621(63) 


1.100 


2.704(60) 


4.8 


1.090 


2.544(61) 


1.105 


2.481(65) 


1.120 


2.623(49) 


4.6 


1.140 


2.549(65) 


1.150 


2.602(64) 






4.4 


1.170 


2.553(51) 


1.160 


2.441(59) 


1.175 


2.612(61) 


4.2 


1.210 


2.548(50) 


1.195 


2.515(62) 






4.0 


1.255 


2.549(60) 


1.225 


2.488(56) 







Table 2: Parameters of (3 and n for fixed v. 
Then we study in Fig. 14 n dependence for some b. One can find that the 
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Figure 14: n dependence of Coulomb term 

action for n > 4 blocked instantons looks on RT. The scaling can be seen more 
clearly in Fig. 15 in which we show the Coulomb term vs. b for fixed v ~ 2.56. 
This is very interesting. 

6. Let us determine the gauge coupling g which is super-renormalizable in three 
dimensions. For that purpose, let us rewrite the action, restoring the lattice 
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Figure 15: Proportional coefficient of Coulomb term 

spacing b = na explicitly: 

S[k] = J2 b6 H bx ) D ( bx -bx')k(bx') 

x,x' 

= Self term + (3 m J2b 6 k{bx) A L \bx - bx')k{bx'), (4.11) 

x,x' 

where f3 m is a constant. Since the mass dimension of the lattice Coulomb 
propagator is 1 and that of S 3 (bx — bx') is 3, the coefficients Cj in (4.7) have 
the mass dimension 2. f3 m is written as 8ir/g 2 . Hence the constant in (4.9) is 
expected to behave as 

Const. ~ ^1. (4.12) 
g l b l 

The 1/b 2 behavior in each gauge is actually observed for b > 1 as shown in 
Fig. 15. ^From the fitting (4.12), we can obtain the renormalized gauge coupling 
g numerically. In the case of v ~ 2.56, we find g = 3.9(2) in the unitary gauge 
condition and g = 3.8(1) in MA gauge. Gauge independence is seen also in 
this case. 

4.5 String tension of the Polyakov model 

Now let us calculate the string tension a c \ = Mg 2 /(An) of Polyakov's model (2.27) 
using the lattice parameters determined from the renormalized trajectory. 
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Note that 

M* = 2^NVW w e-^ X/92) (4.13) 

and 

N=^fa(\/g% (4.14) 

As evaluated in Sec. 2, we get e(oo) = 4ir x 1.67 and cc(oo) ~ 4tt x 3.33. Also we have 
seen above that, for v ~ 2.56, g = 3.9(2) in the unitary gauge and g = 3.8(1) in MA 
gauge. 

Using these parameters, we can estimate the string tension from Eq.(2.27). In 
the case of the unitary gauge condition, 

= 1.4(9) (unitary gauge). (4.15) 

Note that we are considering the scale using a^. The error is not small, since the 
errors of g and v is enlarged by the exponential factor in Eq.(2.16). 
In the case of MA gauge, we obtain 

J— = 1.3(8)(MA gauge). (4.16) 

V &L 

Considering the fact that the saddle point approximation used in deriving Eq.(2.27) 
is not perfectly justified for such small g used here, the above agreement is rather sur- 
prising. In conclusion, we may say that the confinement scenario of three-dimensional 
Georgi-Glashow model by Polyakov is thus verified numerically. 




5. Concluding remarks 

We have studied three-dimensional lattice Georgi-Glashow model in the London limit 
A — > oo under the restriction to an intermediate region where the screening appears. 
In Ref.[14], the screening is seen around the range of 0(2m^/ff). In our case this 
value is around 20 in our unit. As seen in Fig. 15, we are dealing with the region up 
to b ~ 2 in our unit. Hence the screening is not viable in this region. 

We have found that Polyakov's quasi-classical analyses on the basis of the dilute 
instanton gas approximation are in almost quantitative agreement with the Monte- 
Carlo results using the DeGrand-Toussaint definition of a lattice instanton both in 
the unitary and the MA gauges. 

There are other exciting topics to which one can apply the method developed 
here. 

1. When we integrate out non-zero modes of gauge and quark fields perturbatively 
in high temperature QCD in four dimensions, we get an effective model de- 
scribed in terms of the zero-mode bosonic fields alone. The model is equivalent 
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to the three-dimensional Georgi-Glashow model [18]. The effective couplings 
g, v and A are calculated in terms of the quantities in four-dimensional QCD. 
In this case we need to perform simulations for the case with finite A. It is very 
interesting to know if non-perturbative effects in the deconfinement phase of 
(QCD) 4 can be explained totally by the instanton effects. 

2. It is expected that, when A = 0, the three-dimensional Georgi-Glashow model 
can not be described in terms of the Coulomb gas of instantons and the de- 
confinement phase will appear [19]. There is an attractive interaction between 
instantons and anti-instantons which does not depend on the sign of the mag- 
netic charge. It is a Yukawa- type interaction for A ^ 0. But when A = 0, the 
Higgs field becomes massless and the interaction just cancel the Coulomb repul- 
sive interaction between the same-sign instantons or between anti-instantons. 
However, such a scenario needs information of the core region of instantons 
and so non-perturbative numerical evaluation is highly desirable. 
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A. Confinement from the instanton dilute gas model 

String tension from the dilute gas model Eq.(2.27) is calculated as follows. 
Let us start from Eq.(2.24): 

Z[ V ] = J Dxexpj-i {£) J d 3 x[(^(x-^)) 2 -2M 2 cosx]|, (A.l) 

ri(x) = [ dS- (A.2) 
Js \x — y\ 6 

The path integral may be evaluated by the saddle point approximation when the 
coupling g is large enough. Then we have to solve the saddle point equation 

«9 2 (Xci-^) = M 2 sin X ci. (A.3) 

We assume a Wilson loop that the contour C is planar and is placed in the xy 
plane as in Fig. 16 . Then Eq.(A.3) takes the form 
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Figure 16: Wilson loop 

d\ cl = 2ir5(z)6s(xy) + M 2 sin Xcl , (A.4) 

where 

Far from the boundaries of the contour, Eq.(A.4) is approximately one dimensional 
(Xd depends only on z) and has a solution 



i 4arctan(e _Mz ) z > 
Xcl=<! -4arctanfe M n ^<0 ' ( } 



ry(x) only depends on z approximately, too. Eq.(2.25) is rewritten as 

/'' i l X 3 
dyi / dy2 Tt — 2T372 
J-t {{x 1 -y 1 ) 2 + {x 2 -y 2 y + x 3 2 } d / 2 

= arctan 



+ arctan 



+ arctan 



+ arctan 



(r -x 1 )(t-x 2 ) 

X3 \l (x~i - r) 2 + (x~2 - t) 2 + x 3 2 1 
(r + Xi)(t - x 2 ) 
x 3 \J {xi + r) 2 + (x 2 - t) 2 + x 3 2 t 

(r - + x 2 ) 
x 3 ^J ( Xl - r) 2 + (x 2 + t) 2 + x 3 2 

(r + Xl )(t + x 2 ) 
x 3 ^(x 1 +r) 2 + (x 2 + t) 2 + x 3 2 



Now, if r xi and t ^> x 2 ( far from the boundaries of the contour), it is reduced 
to 

n{z)~ 4 arctan (-^^J. (A.7) 
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Hence, one may estimate the integral of Z[rf\ in terms of a z integral alone(where 
x 3 -> z ): 

Zfo] ~ exp |i ^ 2 SM 2 J dz[(xM ~ r]{z)) sin Xcl + cos XcI ]} , (A.8) 

where 5 is the area of the Wilson loop (S = R ■ T = Art). 
/From Eq.(2.25), (A.8) and (A.7), 



\V[C] - oxp^ -^l^:) m2 



dz 



16 arctan 



-Mzi 



1 -e 



-1Mz\ 



:Vr 2 + t 2 + z 2 (l + e -2M,)2 



The string tension is given by 



arctan 



M 2 rt 



e- z (l-e 



-2z\ 



:^{Mr) 2 + (Mt) 2 + 



1 + e- 



-2z^2 



(A.9) 



(A.10) 



The z integral of Eq.(A.lO) is finite if Mr,Mz — > oo . Thus we obtain the string 
tension Eq.(2.27) 

Mg 2 



47T 



(A.ll) 



B. Swendsen's method 

We determine a lattice instanton action from instanton configurations on three- 
dimensional lattice. A theory of instanton is given in general by the following parti- 
tion function 



+oo 



y- III E i < 

x k(x)=—oo j 



-S[k] 



(B.l) 



where S[k] is an instanton action describing the theory. Consider a set of all in- 
dependent operators which are summed up over the whole lattice. We denote each 
operator as Si[k]. Then, the action may be written as a linear combination of these 
operators, 

S[k] = Y / G t S l [kl (B.2) 

i 

where Gi are renormalized coupling constants. 

The expectation value of operator 0[k] is estimated as 



<o[k]>=± m e 

^ \ x k(x)=~oo / 



0[k] e 



-S[k] 



(B.3) 
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Here, we focus on a certain site x' . We define S[k] which contains k(x f ) and rewrite 
the numerator of Eq.(B.3) as 

+00 \ 

II E )0[k\e- s ^ (B.4) 

x k(x)=—oo/ 

+00 \ +00 

II' E \ e-^ G ^- § ^{ J2 0[k] e-^i Gi§i ^}, (B.5) 

x k(x)=—oo/ k(x')=—oo 

where U' means the product except for the site x'. We rewrite Eq.(B.5) as 

(+00 \ 
x k(x)=—oo/ 

+ OO \ 

II E 6[*]e-^°^W, (B.6) 

1 k(x)=—oo/ 

° [k] = v+~ ^Piwi ' (R7) 

fc=— 00 

fc = fc(x') and {A;}' is the coset of k. 

When we consider the DeGrand-Toussaint definition of instanton, the sum with 
respect to instanton charge in Eq.(B.7) is not [—00, +00] but to [— (3n 2 — 1), 3n 2 — 1], 
where n is a step of block spin transformation of instanton. 

Hence we get an identity 

< 0[k\ >=< 6[k\ > . (B.8) 

Since we don't know the correct set of coupling constants Gi, we have to start 
from a trial set of coupling constants Gi . We define 0[k] in which the true coupling 
constants Gi in Eq.(B.7) are replaced by the trial ones Gi as 



where 



mm 



If the trial Gj are not equal to the true G i: < 0[k\ — 0[k] >^ 0. When {Gi} are not 
far from {Gi}, we get 

< 0[k] - 0[k] >~*£<-DSi- (OSi) > (Gi - Gi). (B.10) 
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Actually, if we use Si[k] as an operator 0[k], we get a very good convergence . Hence 
we use here 

< Si[k] - S t [k] >~J2<$i Sj - (S^~) > (Gj - Gj). (B.ll) 

3 

Since the value of Si[k] is obtained from vacuum configurations, we can solve Eq.(B.ll) 
iteratively to obtain the real set of coupling constants Gj. 
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